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To explore conditions underlying the superconductivity in electron-doped TiNCl where T c = 16 
K, we calculate the electronic structure, Wannier functions and spin and charge susceptibilities 
using first-principles density functional theory. TiNCl is the first high-temperature superconductor 
discovered in the a-structure of the layered transition-metal nitride family MNC1 (M=Ti, Zr, Hf). 
We construct a tight-binding model based on Wannier functions derived from the band structure, and 
consider explicit electronic interactions in a multi-band Hubbard Hamiltonian, where the interactions 
are treated with the random phase approximation (RPA) to calculate spin and charge susceptibility. 
The results show that, consistent with TiNCl being a nonmagnetic material, spin fluctuations do 
not dominate over charge fluctuations and both may have comparable impact on the properties of 
the doped system. 



I. INTRODUCTION AND BACKGROUND 

High-temperature superconductivity has been a widely 
pursued subject for condensed matter physicists for 
over two decades. There are several classes of mate- 
rials where unconventional superconductivity is found: 
cuprates whose highest T c 's remain unrivaled, the re- 
cently discovered and intensively studied iron pnictides, 
Nai_ :E Co02 intercalated with H2O, BaBiOs doped by K, 
the Pu-based "115" heavy fermion series, and the tran- 
sition metal nitride halide MNX (M=Zr, Hf; X=C1, Br, 
I). MNX crystallizes in two structures, labeled a and /3. 
The Zr and Hf members were found to be superconduct- 
ing with unprecedented critical temperatures for nitrides 
(15K, 25K) in the /3-structure, which is isostructural to 
SmSI, and contains double-honeycomb layers of alternat- 
ing M and N atoms^. The sister compound TiNCl with a- 
structure has now been discovered to superconduct (16K) 
as wellP 

Based on many examples now, layered structures seem 
to favor high-temperature superconductivity. The re- 
duced dimensionality promotes various instabilities as- 
sociated with Fermi surface nesting. Due to electron- 
electron interactions, many parent compounds of HTSCs 
exhibit long-range magnetic order, such as antiferromag- 
netism (AFM) in cuprates and iron pnictides. The AFM 
order needs to be destroyed upon doping, by electrons 
or by holes, to open the way for superconductivity. This 
is not the case for non-magnetic MNX, which are band 
insulators with a band gap of 2 — 4 eV, and the transi- 
tion metal d states make up most of the lower conduction 
bandeffl 

Both the a- and /3-polymorphs are quasi-2D structures 
with large interlayer spacing and weak van der Waals 
coupling between layers. When doped with electrons, 
A^MNX (A being alkali metals) remains insulating at 
low concentration, then suddenly become superconduct- 
ing at about x = 0.13 in HfNCl and x = 0.06 in ZrNCl, 
and maintain a relatively constant T c up to x = 0.5 [5]. 
The superconducting transition temperature can be as 
high as 26 K, discovered in Li^THF^HfNCP. It has 
also been found that T c may be correlated with the in- 



terlayer spacing, which can be tuned by intercalation of 
different sized molecules 10 . Also, the electron doping can 
be substituted by ion vacancy of the CI atoms with sim- 
ilar superconductivity being founcP^. In a theoretical 
treatment by Bill et al, the dynamical screening ofelec- 
tronic interactions in these materials was modeled 13 * 14 * by 
conducting sheets spaced by dielectrics.^ A fully open 
large superconducting gap without nodes was observed 
with tunneling spectroscopy! 24 " 26 * 

Experimental observations from several perspectives 
confirm that MNX are not electron-phonon BCS super- 
conductors: (1) measured isotope effects are smallp^ (2) 
specific heat measurements^ indicate a small mass en- 
hancement factor; (3) the density of states at the Fermi 
level is small (when electron doped), and T c is almost in- 
dependent of the doping level in the range 0.15 < x < 0.5. 
As a feature peculiar to this system, T c actually in- 
creases as the metal-insulator transition at x cr =0.06 is 
approached, rather than following the common dome 
shape with doping. Linear response calculations also 
agree on the small electron-phonon coupling constant 7 
that cannot account for the observed T c . The impres- 
sive high transition temperatures and easy tunability of 
carrier concentrations (and sometimes effective dimen- 
sionality) suggest there is potential to reach higher T c in 
this class 

The possible candidates of pairing mechanism respon- 
sible for the observed high T c are spin and charge fluc- 
tuations, which have been discussed by both experimen- 
talists and theorists, but opinions remain controversial. 
The specific heat measurement on Li^ ZrNCl is suggestive 
of relatively strong coupling superconductivity, based on 
the observed large gap ratio and specific jump 11 at T c . 
The inter-layer spacing dependence of T c reveals the close 
relation between the pairing interaction and topology of 
the Fermi surface 10 , implying a spin and/or charge in- 
duced superconductivity. The magnetic susceptibility 
measurement of heavily doped LL^THF)^ HfNCl indi- 
cates low carrier density and negligible mass enhance- 
ment factor, in favor of charge fluctuations over the spin 
fluctuations.^ 

A detailed measurement of the doping dependence of 
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specific heat and magnetic susceptibility has been per- 
formed on Li^ZrNCl and the data were compared with 
calculations based on a model HamiltoniarP^, which 
shows correlation between T c and magnetic susceptibil- 
ity. Some change is occurring that affects both, but a 
causal relationship has not been established. There have 
been several band structure calculations for ZrNCl and 
HfNCl that provide the basis for more sp ecific stu dies, 
mostly in the superconducting /3-structure! 1 | 31 | S2 ^ 

Recently, Yamanaka et aP reported superconductiv- 
ity in the alkali metal intercalated a-TiNCl, with T c 
up to 16.5 K. This is the first MNX compound found 
to be superconducting in the a-structure. The lack of 
Fermi surface nesting in a-Li^TiNCl seems to exclude 
large magnetic fluctuations, and its low carrier den- 
sity character is indicative of charge fluctuations as a 
more likely driving force. Theoretically, charge fluctu- 
ation induced superconductivity has been discussed in 
the Hubbard modeP^ES the d- v mo deP^, and been ap- 
plied to study Na^Co02 • yH^CM^iil and organic molec- 
ular superconductors^. A more detailed investigation of 
the electronic structure and dynamical spin/charge sus- 
ceptibility is needed to study the superconducting mech- 
anism in a-TiNCl. In the present work, we conduct a the- 
oretical calculation of spin and charge susceptibilities us- 
ing a many-body Hamiltonian, based on a realistic band 
structure calculated by density functional theory. 



II. CRYSTAL STRUCTURE 

The a-structure^ of the MNX class of compounds, of- 
ten called the FeOCl structure, is shown in Fig. [T] with 
structural data given in Table [j] The a-structure TiNCl 
belongs to space group Pmmn (#59), with 6 atoms per 
unit cell occupying the following sites: Ti(26) (0,|,2Ti), 
N(2a) (|, \,zn) and Cl(2a) (0,0, zci)- The generators of 
Pmmn are two simple reflections x — >• — x and y —y, 
and the non-symmorphic reflection z —> — z followed by 
a (t>, tj,0) translation. 

The Ti-N net within TiNCl is topologically equivalent 
to that of a single NaCl layer. There is strong buckling 
this Ti-N net perpendicular to the b direction, such that 
neighboring chains which are directed along a differ in 
height. These chains are themselves somewhat buckled, 
all of this leading to the placement of Ti ions ±0.8 A from 
the average height, and N ions ±0.4 A from the average 
height. The Ti ions are two- fold coordinated by CI ions 
lying in the y — z plane; the breaking of square symmetry 
of the TiN layer by its strong buckling, can be regarded 
as "due to" this positioning of the CI ions. 

Finally, each Ti is six-fold coordinated by four N and 
two CI atoms. The two Ti-N bonds have very close 
lengths of 2.008 A and 2.015 A, respectively, though the 
N ions lie at different heights in the x and y directions. 
Very roughly, the Ti ion is in octahedral coordination (see 
Fig. 5 of Ref. [8 ), with approximate axes (1,0,0) (toward 
two neighboring N ions), and (0,1,1) and (0,1,-1) (each 



toward one N and one CI ion), and indeed a rough t2 g — e g 
splitting of the Ti 3d states results. The Ti-Ti distance 
is 3.003 A, not much larger than that of Ti-N due to the 
buckled layer structure, so in the tight-binding model we 
construct in the next section, the hoppings between Ti 
sites are also important. 




FIG. 1. Crystal structure of a phase TiNCl (Pmmn, #59). 
Viewpoint is along the b-axis, perpendicular to the buckling 
of the Ti-N net. The buckled Ti-N layer leaves each Ti coor- 
dinated (roughly octahedrally) by 4 N and 2 CI atoms. The 
CI coordinates with Ti along the b-axis, which accounts for 
its orthorhombic structure. 



The experimental lattice constants and atomic posi- 
tions, and relaxed structure parameters with respect to 
total energy, which are used in our calculation, are listed 
in Table [l[ The calculations (see below) confirm the ex- 
pected formal valences. The calculated lattice constants 
are 1-1.5% smaller than the experimental values, but this 
has little effect on the electronic structure. TiNCl is still 
calculated to be an ionic insulator and its theoretical gap 
is similar to what is calculated using the experimental 
lattice parameters. 





a 


b 


c 


ZTi 


ZN 


zci 


Expt. 
Theory 


3.938 
3.891 


3.258 
3.200 


7.800 
7.699 


0.1011 
0.1011 


0.0509 
0.0522 


0.3322 
0.3384 



TABLE I. Lattice constants (in units of A) and internal struc- 
tural parameter z for the three atoms. Experimental values 
are from Ref. 8 The theoretical values are our optimized 
values. 



3 



III. BAND STRUCTURE AND WANNIER 
FUNCTIONS 

A. Methods 

The band structure has been computed by the full- 
potential local orbital minimal basis set method im- 
plemented in the FPLO code!^ The exchange correla- 
tion is treated by the generalized gradient approximation 
GGA96p^and the k-mesh used is 16 x 16 x 8. The effect 
of spin-orbit coupling is small so calculations were done 
in the scalar relativistic scheme. 
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FIG. 2. (color online) Band structure of TiNCl along 
the orthorhombic symmetry lines, calculated with the GGA 
exchange- correlation functional. The thick "fatbands" are the 
tight-binding representation determined by the Wannier func- 
tions labeled in the inset, chosen to represent accurately the 
bands at and around the Fermi level after doping. 



B. Electronic Structure 

The calculated band structure of pristine TiNCl is 
shown in Fig. [2] and is generally consistent with that 
presented by Yamanaka et alP* plotted along other lines 
in the zone. It is an insulator with a calculated energy 
gap of 0.5 eV. The real band gap may be as large as 1 eV, 
based on the common observation that LDA and GGA 
underestimates gaps in insulators. The band structure 
exhibits clearly a two-dimensional feature, gauged from 
the general flatness of bands along the T — Z direction 
perpendicular to the layers. The states on either side 
of the gap are very two-dimensional, considering the ex- 
treme flatness of those bands along T-Z. 

The twelve valence bands are made of six N 2p and six 
CI 3p states, and the conduction band is comprised of ten 
Ti 3d states. The 3d bands show a "t2 g — e g " crystal field 
splitting (three states below and two above), that arises 
in spite of the nonequivalence of the five 3d orbitals in 
this structure. As can be seen in the partial density of 
states plotted in Fig. [3j there is 3d weight in the valence 




Energy (eV) 

FIG. 3. (color online) Total and partial density of states of 
TiNCl, showing the Ti 3d - N 2p mixing. The inset figure is 
a blown-up of the region near Fermi level. 



bands and N weight in the conduction bands, reflecting 
substantial N 2p - Ti 3d hybridization in addition to the 
ionic character reflected in their formal charges. 

Whereas the 3d environment appears locally to be 
pseudo-cubic, the low site symmetry severely splits the N 
2p states, with 2p x and 2p y becoming quite distinct. The 
top valence band is primarily N 2p x character, which ex- 
tends down to —5 eV. The N p y and p z bands have their 
maximum 1 eV lower, and the CI 2p weight is concen- 
trated at the bottom of the valence bands. 

The inset in Fig. [3] shows an enlargement of the to- 
tal and atom-projected DOS around the Fermi energy. 
The onset at 0.5 eV and the smooth slope to 1.2 eV is 
characteristic of a two-dimensional band which becomes 
non-parabolic away from the band edge, and strongly so 
in the 1.2-1.5 eV region. At 1.5 eV the onset of the sec- 
ond band, with its much heavier mass, is clear. However, 
the DOS does not have the sharp step at the top of the 
valence band that is characteristic of a 2D system. 

The Ti 3d orbitals are lifted in degeneracy entirely by 
the orthorhombic point group site symmetry, but as men- 
tioned above the conduction bands are separated by a 
crystal field analogous to cubic t<i g — e g splitting. Check- 
ing the band character reveals that, in terms of 3d or- 
bitals expressed in terms of the orthorhombic coordinate 
axes, d xy , d z i^d xz have most of the weight in the 0.5 — 3 
eV region, and d yzi d x i_ y i are in a higher energy window 
of 4 — 6 eV. Thus it is feasible, in a low-energy tight- 
binding model, to include only Ti d xyi d z i,d xz and N 2p x 
states. Plotted on top of the DFT bands in Fig. [2] is 
the tight-binding fit using the Wannier functions. The 
representation of the full t^g complex is excellent, as is 
that of the top of the upper valence band. 

The distance between TiNCl slabs and the weak inter- 
layer coupling allows intercalation of alkali atoms, which 
act as electron donors. This feature validates the rigid 
band shift approximation in simulating doping. Doped-in 
electronic carriers will go into the single lowest conduc- 
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tion band, which is quite two-dimensional as mentioned 
above but is dispersive within the plane. This band has 
strong Ti 3d xy character, similar to the in-plane Ad char- 
acter in ZrNCl. The Fermi surface of electron-doped 
TiNCl is an oval centered at the V point. This point 
has some relevance for the superconductivity, since with 
a single Fermi surface there can be no nesting of discon- 
nected Fermi surfaces, such as are proposed 12 to play an 
important role in many other layered superconductors, 
such as Fe pnictides, as well as f3- structured ZrNCl and 
HfNCl. The similar characters of TiNCl and ZrNCl, and 
their similar values of T c , suggest that possible nesting 
of Fermi surfaces is not an important feature for pairing 
in the materials. 



C. Wannier Functions 

Because the susceptibilities we will calculate have a 
number of local orbital matrix elements equal to the 4th 
power of the number of orbitals retained, we have calcu- 
lated selected low-energy Wannier functions (WFs) that 
will be used to construct our many-body Hamiltonian, 
using projections of the Bloch states onto the correspond- 
ing atomic orbitals. The four atomic orbitals mentioned 
above allow us to reproduce the bands on either side of 
the gap: Ti d xy , d z i, d xz and N p x . While the Ti u t2 g " or- 
bitals are not optimal in diagonalizing the local "octahe- 
dral" symmetry, they and their relation to the 2p x orbital 
are more readily visualized. Since the RPA calculations 
described below were performed in the electron-doped re- 
gion where the Fermi level is shifted into the conduction 
bands, considering only the N p x WF in the valence band 
is sufficient to understand the q-dependence. Aside from 
being farther removed in energy, the remainder of the 
valence bands form a complex of bands spread uniformly 
over the zone, contributing little to any q-dependence. 

The Ti-N layer is strongly buckled and there are 2 Ti 
and 2 N atoms per unit cell with different z coordinates, 
the actual tight-binding model contains 8 bands and 8 
WFs, but WFs on symmetry related ions are symmetry 
equivalent. Overall the Wannier orbitals generate a well 
represented band structure compared to the DFT bands 
within the energy window of interest, as shown in Fig. [2] 

The hopping amplitudes of the Wannier orbitals are 
listed in Table [TTJ Hopping integrals smaller than 0.05 
eV were not listed because they only marginally alter 
the band structure and obfuscate interpretation. The 
on-site energies of the "£2#" orbitals are 2.25 ± 0.08eV, 
lying within the largest peak of the DOS. The p x energy 
is — 3.22eV, in the middle of the valence bands. Thus 
there is a 5.5eV separation of valence and conduction 
band centers, and a gap of 0.5eV. 

The dispersion within the pair of valence p x bands 
is represented largely by hopping between neighboring 
p x WFs (recall, the p x WF contains some d character, 
and vice versa), both being about 0.5eV. Hopping am- 
plitudes to the d orbitals are t pd ~ 0.2 — 0.3eV. In the 



conduction bands, the d xy orbital has hopping amplitude 
\t\ ~ 0.17eV to its partner within the cell as well as to 
its replicas in neighboring cells in both directions. The 
hopping to the p x orbital (0.33 eV) is twice as large, 
and apparently is the dominant contributor to the 3.5eV 
bandwidth. Due to the relative orientations, hopping to 
the other d WFs is no more than half as large as the 
d X y — d xy one. The other two d WFs form rather nar- 
row bands, reflected by smaller hopping amplitudes; note 
that both have hopping to the p x orbital of 0.21 — 0.24eV. 





[0,0,0] [1,0, As] [0,1, As] [1,1, As] 


[2,0,0] [0,2,0] 


[dxy 7 dxy ) 


2.33 


-0.16 


-0.17 -0.18 


(d xy ,d z 2) 




-0.08 




{d X y 5 d X z ) 




0.07 


-0.11 


(d z 2,d z 2) 


2.18 


-0.05 




(d z 2 , dxz) 




0.19 


0.07 


{dxz 7 d X z) 


2.25 




-0.05 0.13 


{d X y 5 Px ) 


-0.33 






{d z 2,p x ) 


0.24 0.22 






{dxz ? Px ) 


0.21 






(Px,Px) 


-3.22 


0.51 


0.47 0.14 



TABLE II. Nearest, 2nd, and 3rd neighbor hopping integrals 
in units of eV. The hopping vectors are in units of [a/2, b/2, c], 
and Az represents the difference between the z coordinates of 
the two orbitals. This representation is purely two-dimension 
(no coupling along the c axis). 



IV. MANY-BODY HAMILTONIAN AND 
RANDOM PHASE APPROXIMATION 

The random phase approximation (RPA) applies an 
interaction to the non-interacting Hamiltonian 

H = J2H* b cl Ab , (1) 

k,a6 

where a, b are composite orbital and spin indices of the 
basis Wannier orbitals. The interaction Hamiltonian, in 
general, is the following symmetric form 

H l = \ Y U *bcd4 a C ib c\ c C id + Y^ Y V ^cdc\ a C ih c) c C jdl 

i abed (* 5 i) abed 

(2) 

where U and V represent on-site and inter-site (only 
nearest neighbors) interactions, respectively. Q can be 
Fourier transformed into 

Hl = J^YY F ^d(q)4,a C k+q,6 C p+q,c C prf' ( 3 ) 
kpq abed 

in which F abcd {q) = U abcd + j{q)V abcd is the interaction 
kernel matrix. In the second term, 7(q) = ^ l e zq Rz 
(I running over the nearest neighbor pairs of sites) is 
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the structure factor, which brings in q-dependence into 
Fabcd(q)- The bare susceptibility is calculated as 

Xabcdfau) = ^G ad (k,cj)G c6 (k + q,cj), (4) 

k 

where G a 6(k,u;) is the non-interacting Green's function 
MM) = V^», (5 ) 



n 



U + fJL ~ £ nk 



and the summation is taken over all bands. Applying 
RPA, which sums up the higher order diagrams in the 
geometric series, we have the full susceptibility presented 
in a matrix equation 

X(q,w) = [I + x°(q,w)Ite^(q)]" 1 X (q,w) (6) 

where a matrix x is formed from Xabcd by contracting the 
first pair of indices and the last pair of indices. 

So far we have derived a very general formula for an ar- 
bitrary interaction Hamiltonian. To study our case, next 
consider a more specific model in the form of an extended 
Hubbard Hamiltonian, following Kuroki's modeP^, but 
add an extra inter-site interaction term: 



jn ih(J ' (7) 



-J^2S ia - Sib + J' ^2 C la^4al C ibl C m 
a^b a^b 



in which a, b are orbital indices, z, j are site indices of 
the lattice, and a is the spin index. U is the intra-orbital 
Coulomb repulsion, U' is the inter-orbital Coulomb inter- 
action, is the hopping between Wannier orbitals, V a b 
is the inter-site Coulomb interaction between orbitals a 
and 6, J is the Hund's rule coupling, and J' is referred to 
pair hopping between orbitals. From this Hamiltonian, 
the susceptibility is calculated by 



x 5 (^,q) 



x°(w>q) 



x c (^,q) = 



S(q)xV,q)' 



(8) 



J + C(q)x°(a;,q)' 
The interaction matrices S(q) and C(q) take the form: 

[/, a = b = c = d 

U' , a = c b = d 

J, a = b ^ c = d 

J 7 , a = d b = c 

U + 2Kc Ke 7(q), a = b = c = d 
-U' + J, a = c^b = d 



S a bcd — 



(9) 



Cabcd 



2U' - J + 2V r ac Re7(q), a = 6^c = d 



where [/, t/ 7 , J, and J 7 terms appear only if all indices 
are orbitals on the same site. Finally, we can also cal- 
culate the macroscopic susceptibilties by performing a 
summation over the orbital indexes: 

X moc (q,w) = ^SijXiM^Skl, (io) 

in which 5^ is overlap matrix, and in our case it has the 
form of delta-function 5**. 



V. SPIN AND CHARGE SUSCEPTIBILITY 

With our model just constructed, we calculate the spin 
and charge susceptibilities. The model is a multi-band 
extended Hubbard model on a 2D rectangular lattice 
with 4 sites (two Ti and two N) per unit cell. For the 
on-site interaction terms, we use Udd = U' dd = 1.5 eV, 
U pp = 1.0 eV, and J — J' = 0.2 eV. These values are 
somewhat smaller than might be used in a traditional 
Hubbard model calculation; this is partly to compensate 
for the fact that RPA has a tendency to overestimate the 
strength of the interaction due to the lack of the self- 
energy correction. 30 Moreover, we are using WFs rather 
than atomic orbitals, for which the extension onto neigh- 
boring sites will suppress the intra-atomic interactions U 
and J. 

For inter-site interactions, we assume V ac to be spin 
and orbital independent, and that it only depends on 
the distance between the two sites. Taking into account 
the Ti-N and Ti-Ti distances mentioned above, we use 
Vti-n — 0-5 eV (nearest neighbor), Vti-ti — 0.3 eV 
(2nd nearest neighbor). Since the Wannier functions 
have contributions from neighboring sites, it is reason- 
able to set the inter-site Coulomb repulsion V slightly 
larger than traditionally used for atomic orbitals. The 
calculation is done at T = 0.02 eV (220K) and u = 0, 
with a k-mesh of 40 x 40 x 4 and q-mesh of 20 x 20 x 2. 
The occupation is set at 4.3, simulating x — 0.15 electron 
doping in Ao.isTiNCl (since there are two formula units 
per unit cell) by raising the Fermi level into the lowest 
conduction band. 



O 



J', 



d^b 



FIG. 4. (color online) The spectral function — ImG(k, uj 
in the k x — k y plane, showing the oval (nearly circular) F 
surface for x=0.15 electron doping. 
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Figure [4] shows the magnitude of the imaginary part 
of Green's function, which provides a view of the Fermi 
surface. With a simple, nearly circular Fermi surface like 
this, the bare susceptibility is expected^ to be isotropic 
out to q = 2k p, with a relatively constant plateau be- 
havior inside 2kp radius. The inter-site Coulomb inter- 
action can give rise to charge fluctuation, creating collec- 
tive electron motion and possible charge ordering. Com- 
petition between on-site and inter-site interaction of d 
electrons can lead in principle to a frustration of both 
spin and charge ordering. The hybridization between d 
and p orbitals opens another channel, that of a charge 
transfer instability. One of the interesting questions is 
whether some combination of these processes can create 
excitations that can pair up electrons, analogous to the 
behavior found in a d-p model in the limit of infinite U 
and nearest neighbor hybridization 21 . 




(c) 

FIG. 5. (color online) Representative orbital spin suscepti- 
bilities. In the basal plane, the axis on the left (longer in 
perspective) is q x axis, and the one on the right (shorter) is 
q y axis, (a) xfiii (b) X1133 ( c ) Xim- Note that the maxi- 
mum in X1133 i s located at (7r/a, 0, 0) and not at 2/cf- Note 
the different vertical scales on the panels. 

In Figure [5] some representative spin susceptibilities 
in the orbital representation are plotted on the two- 



dimensional basal plane (q z = 0) in the BZ. To clarify, we 
denote the orbitals by numbers in the order: (l)Tii-d xy , 
(2)Ti 2 -d xy , (3)Ti!-^ 2 , (4)Ti 2 -d*2, (5)Th-d xz , (6)Ti 2 - 
d X zi (7)Ni-p £C , (8)N 2 -#e. The largest spin susceptibility 
is found for X1111 (intersite, xy O xy) which has approx- 
imate 4-fold symmetry for magnetic fluctuations of the 
same orbital d xy on Ti sites. The anisotropic behavior 
of 

X1133 (on-site, xy z 2 ) comes from the orthorhom- 
bic symmetry of the lattice, i.e. a ^ 6, which brings in 
anisotropic q-dependence, in this case strongly so. The 
spin fluctuations between d and p orbitals Xim have a 
sizable overall magnitude, comparable to d — d fluctua- 
tion, but small variation with q, because four neighboring 
N atoms have almost the same distance to Ti. Due to 
the lack of Fermi surface nesting, there is no divergent 
behavior in the spin susceptibility, presenting different 
physics from the /3-HfNCl which has two circular Fermi 
surfaces located at two high-symmetry (K) points in the 
BZ which can provide near perfect nesting. 




(b) 

FIG. 6. (color online) Representative charge susceptibilities 
in the full Brillouin zone, (a) xfiii (b) Xim- The latter 
quantity is hardly distinguishable from its spin counterpart. 
Note the different vertical scales on the panels. 

Representative charge susceptibilities are shown in Fig- 
ure [6] Often they have similiarities to the spin suscep- 
tibilities, with comparable but somewhat smaller magni- 
tudes. It is known that in an extended Hubbard model on 
a square lattice, at zero frequency, x S an d X° have simi- 
lar q-dependence and only vary in magnitude.^ Without 
the long-range Coulomb interaction, charge fluctuation 
will always be smaller than spin fluctuation because of 
the different signs in the RPA formula. The difference 
between x S an d X° wm become more apparent at non- 
zero uj. The q-dependence of x C is similar to that of x S 
but shows somewhat more structure in Xmi- Note that 
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(c) 



FIG. 7. (color online) Macroscopic susceptibility: (a) bare; 
(b) spin; (c) charge. Note the strong deviation from square 
symmetry. 

the magnitude of Xmi 1S om y na ^ that of Xmi- The 
spin and charge fluctuations within the N p-channel are 
very small since the p orbitals are fully occupied so fluc- 
tuations only happen as a second order effect. However, 
the presence of the N p band very close to the lowest 
d conduction band opens an additional channel for fluc- 
tuations between them (xfm? Fig. p5p and Xu77> Fig. 

To close the comparison, we show the macroscopic sus- 
ceptibilities in Figure [7| The imprint of 2k f is evident. 
Beyond 2kp, the susceptibilities decrease slightly and 
with (near) square symmetry. Inside 2kp, the variation 
is greater and displays the rectangular symmetry of the 
lattice. The overall spin enhancement of the macroscopic 
susceptibility (x S (o) / X° (o)) near q = is about 1.4. For 



the model with the realistic parameters considered here, 
our test calculations show that x S wm approach diver- 
gent behavior when U ~ 4 eV. Since TiNCl has wide 3d 
bands, it is unphysical to use U anywhere near 4 eV (we 
are using 1.5 eV), so strong spin fluctuation is unlikely 
to occur in this material. The macroscopic charge sus- 
ceptibility has smaller magnitude and q-dependence, but 
will show divergent behavior when inter-site interaction 
V is much larger than U, but that regime is unrealistic 
for the case studied here. Overall, both spin and charge 
susceptibilities show moderate enhancements compared 
to the bare susceptibility, without any approach to an 
instability toward spin or charge ordering. 

VI. CONCLUSION 

To conclude, we have constructed a many-body ex- 
tended Hubbard model using a realistic band structure 
obtained from density functional theory calculations. 
RPA is applied to obtain the spin and charge suscep- 
tibilities. In a system like a-TiNCl, where the crucial 
ingredients of high-temperature superconductivity, such 
as strong electron-phonon coupling and good Fermi sur- 
face nesting, seem to be missing, spin and charge fluc- 
tuations are the remaining candidates. Our calculations 
show that the spin and charge enhancements of suscep- 
tibilities, both intra-band and inter-band, are small due 
to moderate correlations. However, spin and charge fluc- 
tuations can produce substantial values possibly capable 
of encouraging electrons to pair. Although spin fluctua- 
tions are present in a-TiNCl, it is worth pointing out that 
the physics is very different from the /3-structure coun- 
terparts even though both are nonmagnetic, and appar- 
ently distinct from the recently discovered Fe pnictides 
where magnetism is an important feature in parent com- 
pounds. As the nonmagnetic nature of TiNCl and other 
MNX materials indicates, as well as seen from the results 
from our RPA calculation, charge fluctuations may have 
an important role in superconductivity in these systems. 
Overall, we still do not have a clear understanding of how 
superconductivity arises from the fluctuations, as with all 
other high-Tc families. 
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